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Abstract 

We prove a new global existence result for the asymptotically flat, 
spherically symmetric Einstein- Vlasov system which describes in the 
framework of general relativity an ensemble of particles which interact 
by gravity. The data are such that initially all the particles are moving 
radially outward and that this property can be bootstrapped. The 
resulting non-vacuum spacetime is future geodesically complete. 
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1 Introduction 



Global properties of non-vacuum spacetimes are one of the central themes 
of the mathematical analysis of general relativity. A matter model which 
has proven particularly suitable for such an analysis is the collisionless gas 
as described by the Vlasov equation. Here the matter content of spacetime 
is represented by a large ensemble of particles with a density function / on 
phase space, i.e., on the mass shell of the tangent bundle of the spacetime 
manifold. The particles move along time-like geodesies which is equivalent to 
/ satisfying the corresponding continuity equation, i.e., the Vlasov equation, 
and the energy momentum tensor is then expressed in terms of /. 

For a general introduction to the resulting Einstein- Vlasov system we 
refer to [1] and [TO]. A major result for the asymptotically flat case of this 
system is a global existence result for small, spherically symmetric initial 
data [5]. In addition, there is numerical evidence [2j [5j [9] that at least in 
the spherically symmetric case the weak cosmic censorship conjecture holds: 
The gravitational collapse of regular initial data leads to a black hole rather 
than a naked singularity. 

In view of the small data result mentioned above and the lack of proof of 
global existence for general data (in Schwarzschild coordinates or in maximal- 
areal coordinates) it is natural to ask if one can find an open set of initial data 
which are not small in the sense of [5] but for which global existence holds. 
In the present paper we do find initial data with this property by establishing 
global existence for data which are such that initially all the particles in the 
ensemble move radially outward sufficiently fast so that they can be shown 
to keep on moving outward for all future time. This configuration results 
in estimates on the geometry of spacetime which imply that it is future 
geodesically complete. 

The basic mechanism of the proof is completely different from the one 
for small initial data which relied on the dispersive properties of the Vlasov 
matter model when the fields are small. Since its actual implementation for 
the Einstein- Vlasov system quite hides this mechanism, it should be useful to 
explain it first in the case of the much simpler Vlasov- Poisson system which 
is the Newtonian limit of the Einstein- Vlasov system. Due to the spherical 
symmetry the maximal force experienced by a particle at distance r from the 
origin is —M/r 2 in the Vlasov-Poisson case, where M > is the total mass 
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of the ensemble. Hence along any particle trajectory 



d (\ 2 M\ . , M . ( . , M 



as long as its radial velocity r = ty = x- f/r>0. Hence 

1 9 . , M 1 , M 
V(t) — > -w 2 (0) 



2 v ' r(t) ~ 2 w r(0) 
and 

1 „. , 1 9 . . M 
-w 2 (t) > -w 2 (0) 



2 w ~ 2 v ' r(0) 

on any time interval on which w(t) remains non-negative. Now let u>_ > be 
an initial lower bound for the radial velocities of the particles in the ensemble, 
r_ > an initial lower bound for their distance from the origin, and assume 
that 

1 9 M 

W_ := -w 2 > 0. 

2 ~ r_ 

Then as lone; particle is moving outward, 

w(t) > W-, r(t) > r_ + W-t. 

But this implies that all the particles keep moving outward for all future time. 
For the Vlasov-Poisson system global existence and uniqueness of smooth 
solutions has been established for general initial data [BIS], but the above 
argument can easily be extended to yield global existence for spherically 
symmetric initial data as specified above. 

For the spherically symmetric Einstein- Vlasov system we find a substitute 
for the pseudo-energy |u> 2 — — used above, but when taking its derivative 
along a particle trajectory it turns out to be important to use suitable co- 
ordinates on spacetime, namely maximal-areal coordinates, in order to see 
that under suitable circumstances this quantity does not decrease along radi- 
ally outgoing trajectories. Once this is shown a bootstrap argument implies 
that all the particles keep moving radially outward on the future directed 
maximal interval of existence of the solution. We can then use a result from 
[TO] to conclude that the solution is global in the time coordinate used, and 
the obtained control of the metric quantities turns out to be sufficient to 
establish future geodesic completeness. 
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2 The Einstein- Vlasov system 

We consider the system in the spherically symmetric, asymptotically flat case 
and write the metric in the following form: 

ds 2 = -(a 2 - a 2 (3 2 )dt 2 + 2a 2 (3dtdr + a 2 dr 2 + r 2 (d6 2 + sin 2 9 d(p 2 ) . 

Here the metric coefficients a, (3, and a depend on t e R and r > 0, a and a 
are positive, and the polar angles 9 G [0, tt] and e [0, 27r] parameterize the 
unit sphere. The radial coordinate r is thus the area radius. Let K% be the 
second fundamental form and define 

K = Kg = A. (2.1) 

By imposing the maximal gauge condition, which means that each hyper- 
surface of constant t has vanishing mean curvature, we obtain the following 
field equations where subscripts denote partial derivatives with respect to r 
or t: 

a r = ^a 3 rK 2 + 4irra 3 p+ ^-(1 - a 2 ), (2.2) 
3 

k t = — k — Attclj, (2-3) 
r 

a t = 2aan + (aj3) r , (2.4) 

(?L - l\ + ^ ( 2r ^ + a 2 - l) + 47ra 2 a(S - 3p). (2.5) 
\a r J r l \ a J 

The Vlasov equation takes the form 

« + d - ! + " *»» + £e) a ~ s = °< < 2 - 6 > 

where 



E = E(r, w,L) = a/1 + u> 2 + L/r 2 . (2.7) 

The variables w and L can be thought of as the momentum in the radial 
direction and the square of the angular momentum respectively, see [7j for 
more details. The matter quantities which appear in the field equations are 
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defined by 



oo />oo 



p(t,r) = ^ I I Ef(t,r,w,L)dLdw, (2.8) 



oo JO 

OO POO 



j{t,r) = - 2 I I wf(t,r,w,L)dLdw, (2.9) 

T J_ 00 J 

/oo /*oo 2 i t I 2 
/ ^ + F /r f(t,r,w,L)dLdw. (2.10) 
■oo JO ^ 

We impose the following boundary conditions which ensure asymptotic flat- 
ness and a regular center: 

a(t, 0) = a(t, oo) = a(t, oo) = 1. (2.11) 

The equations fl2.2p - fl2.10p together with the relation (12. ip and the bound- 
ary conditions (12.111) constitute the Einstein- Vlasov system for a spherically 
symmetric, asymptotically flat spacetime in maximal-areal coordinates. 

First we note that the phase space density / is constant along solutions 
of the characteristic system 

, = a^>_ 
a(r, r)h 

a r (r,r)E , . . . a(r,r)L 

™ = ~ f \ ~ Mr,r)K(T,r)w+ / 2.13 
a(r, r) a{T,r)r 6 E 

L = 0, (2.14) 

of the Vlasov equation. If r i — ^ (i?, VF, £)(t, t, r, to, L) denotes the solution of 
the characteristic system with (R, W, L)(t, t, r, w, L) = (r, w, L) then 

f(t,r,w,L) = f((R,W,L)(0,t,r,w,L)), 

o 

with / = f\t=o the initial data for /. In particular, / > provided this is 
true for the initial data as we will assume throughout. Next we notice that 
(12.31) can be rewritten as 

{r 3 Kj r = —4:7ir 3 aj, 

and upon integration 

47r r 

K (t,r) = t a(t,s)j(t,s)s s ds. (2.15) 
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The Hawking mass m is given by 



We also introduce the quantity 



(2-17) 



and note that by (12.21) . p can be expressed as 



p(t,r) 



J \ 4np(t,8) + ^K 2 (t,s)^ s 2 ds. (2.18) 

Assuming that the matter is compactly supported initially and hence also 
for later times, (I2.15P implies that n(t, r) ~ r -3 for r large. Hence the limits 
as r tends to oo of m and p are equal, so that the ADM mass M is 

M = J ^4vrp(t, r) + ^ 2 (t, r)^j r 2 dr. (2.19) 

The ADM mass is conserved, and p(t,r) < M. 

Finally we notice that using (12.21) the second order equation (12. 5p can be 
rewritten in the form 



— a r = 47rr aa(p + S) + 6r aan , 
a 



and upon integration, 

aft r) f r 

a r (t,r) = 2 / (4vraa(p + S) + 6aa/t 2 ) s 2 ds. (2.20) 

Thus a is monotonically increasing outwards, and from (12.111) it follows that 

< a < 1. (2.21) 
A direct computation using (|2T5]1 . (P22D, fl2TT7D . and (12TT81 implies that 

\ ^2r^ ^22 yiS ^32 

—a r — ap = Attt ba H — r k a — Anr pa r r k a r . 

a 2 J r 2 2 



With the boundary conditions (a/i)| r=0 = = (r 2 a 2 av)| r . =0 this yields the 
identity 



a r (t, r) = (afi)(t, r) + 



9 3 
ArrSa + -^k 2 ol — Affspa r — -SK 2 a r 



s 2 ds, 
(2.22) 



a 2 (t,rj Jo 
which turns out to be very convenient in what follows. 

3 The main result 

o 

Let / be some initial data of ADM mass M such that 

o 

supp/c]r_,r + [x[w_,oo[x[0,L + ] . 

We wish to specify these parameters in such a way that for the corresponding 
solution of the Einstein- Vlasov system all particles have radial velocity w > 
with an additional parameter W_ > 0. To this end we abbreviate 

£ = — = = (3 1) 

r- ' V y/T+WT+Ljrl' ' 

and we require that the above parameters satisfy the following conditions: 



2 mm { — , , 77 > — 1 , 7T W 



and 



— (Wf + ^l + ^jj. (3.3) 

After stating our main result we will show that there exist parameters satis- 
fying these conditions, and that there exist initial data with such parameters. 

Theorem Suppose that < r_ < r + , W- > 0, L + > 0, W- > 0, and 
M > satisfy the conditions A3.2\) and A3.3\) , and consider a solution f on 

o 

a maximal existence interval [0,T[ such that the initial data f = f\ t =o is 
compactly supported, smooth, non-negative, has ADM mass M, and satisfies 

supp/c]r_,r + [x[w_,oo[x[0,L + ] . (3.4) 



7 



Then T = oo, and the resulting spacetime is future geodesically complete. 
The proof of this theorem is given in Section H] below. 

To show that parameters £, r ± , w_, L + , W_ exist which satisfy (13.21) and 
(13.30 we proceed as follows. Let 

1 — .V Q P 3 P 

Since D(0, 1) = 1 and D(l/3, 1) < there exists some £* G]0, l/3[ such that 
1) > for £ G [0,f*[ and £>(£*, 1) = 0; numerically, « 0.0994. So 
firstly we assume that 

Secondly, we choose r + > r_ > such that 

r_ 9 £ 3 £ 



where we notice that the right hand side is increasing in £ and its value for 
£ = £* is ~ 1.405 so that the inequality can be satisfied with r + > r_ > 0. 
Now we fix L + > arbitrarily. Since rj — > 1 as W 7 - — > oo we can choose 
W_ > sufficiently large so that 



and 



o r - , 9 £ 3 £ 
2— > 1 + - + - -7== + 



>+ ' ' 4 l-£ ' 2 v 7 !^ ' ' 

Then the factor multiplying W£ in the condition (13.21) is positive, and by 
further increasing W_ if necessary condition (13. 2p is satisfied. Finally, we 
take u>_ > sufficiently large so that (13. 3p holds as well. 

Let us now fix a set of parameters £, r±, L + , u>_, W_ such that (13. 2\ and 
are satisfied. We want to argue that an initial condition with these pa- 



rameter values does exist. In order to do so we fix some compactly supported, 
smooth function g such that 



supp g C]r_, r + [x oo[x [0, L 



+J > 



and we consider /= v4g for some amplitude A > 0. Quantities induced by / 
will be denoted by the corresponding superscript, and they depend on A. By 
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continuity, we can make sure that d\r) < 2 for all A > sufficiently small; the 
dependence of the initial data for the metric quantities on the initial matter 

o 

terms and hence on / will be studied in |3] . Now let 

POO 

1= 4nfcs)s 2 ds. 
Jo 

Then by f l2TT5|) . 

I*r)|<^ [ r d{s)\j(s)\s 3 ds<% [\s)s 2 ds = ^. 
Since «(r) = for r < r_ this implies that 

f°° ( 3 \ fif 2 

I<M = J Uirflr) + ^\r) J r 2 dr < I + 

o 

Since the quantity / depends linearly on the amplitude A, 2M(/)/r_ = ^ < 
^* for A > sufficiently small, and arguing again by continuity any value of 

o 

£ e]0,£*[ should be realized by adjusting A properly. Since / has the same 

o 

support as g, all conditions required in the theorem then hold for /. 



4 Proof of the theorem 

Step 1: The bootstrap argument. 

Our first and crucial step in the proof of the theorem is a bootstrap argument 
which shows that the condition w > W_ which holds on the support of / 
initially will persist on the maximal forward existence interval of the solution. 
To begin with we show that for all t G [0, T[ and r > r_, 

1 1 — 3£ 
a(t,r) < a + = and a(t,r) > a_ = -; (4.1) 

v 1 — 4 1 — s 

eventually we shall see that these bounds hold for all r > 0. We recall that 
£ = 2M/r_. By (EUD, 

1 1 2»{t,r) >1 2M 



a 2 (t, r) 
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for r > r_; the inequality is strict because the solution is non-trivial and 
hence re does not vanish identically. By (I2T2TID . (I2T2T|) . and (gUgj) , 



a. 



•(*,0<-±/ {^{p + S) + QK 2 )s 2 ds<^-. (4.2) 

Since a(t, oo) = 1, this implies that 

, , , 4a2M l-3£ 
a(t,r)>l--±- = T -^. 

Next we define for t G [0, oof, 

R - * = + — 71 Jr r , r (4-3) 

a+ y/1 + W 2 + L + /r 2 _ 
iMf)=r + +(l + ^)t, (4-4) 

as well as 

. fl_(t) . fr_ a_ W_ 

7 = mm — = mm 



e[0MR+(t) [ r+ a+(l + a+M/rJ) + Wl + L + /r 2 _ 

. / r l-3£ \ 
mm < — , — — — 77 > , (4.5) 



r + 'v^e + e/2 

where we recall the abbreviations in (13. ip and (14. ip . Let [0,t [ denote the 
maximal time interval such that for t G [0, to[ and (r, w, L) G supp fit)-, 

R„(t) < r < R + (t), w > W_ . (4.6) 

By (13.31) and (13.41) this holds for t — 0, and by continuity there is a maximal 
time interval [0,t [ with to > on which (14. 6p persists. Suppose that t G 
]0, T[. Then we must have equality at t = to in at least one of the inequalities 
in ( 14. 6p . But in the following part of the proof we show that (14.61) persists 
on [0, to]- Thus necessarily t = T. 

Consider a characteristic such that f(t,r(t),w(t),L) > 0. Denoting the 
parameter along the curve by t it is calculated that 



d_ 

dt 



2M 



E 2 



, , , 2M\ a r _, B 2 M 2 

-2 1 -£w - 4-w 2 + 6—Bw 2 

r J a r r 1 

L M , Ma .1/ , , „. 

+ 2—B — 6—BL + 2——Ew — 2— (3, 4.7 

iyO f* f* Oil f* 
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and we want to show that this quantity remains positive on the time interval 
[0, t„[. By (J21D and (Oil , j > for i e [0, t [- Hence (EH) and (T^T^j) lead to 
the estimate 



/3(t,r) 



a(t, r) 



/■r 

/ 47r a(i, s)s 3 ds < 0, 
Jo 



(4.8) 



and we can drop the fifth and seventh term in (14. 7p . Using the estimate 

1 _ 1 2/i(t,r) >i 2M 
a 2 (t,r) r ~ r 

and the fact that a r . > 0, cf. (I2.20p . we find that 

2M 



d 
dt 



E 



> -2-^Ew - A-w 2 + -rr I 3Mw 2 + - I + 2—-Ew. 
a 6 r r 2 \ r / r 2 a 

(4.9) 



By (I2~22j) and in view of a(t, s) < a(t, r), 

M. 9 . . , (aa 2 )(t,r) 
«r(t,r) < ^(aa 2 )(t,r) + ± -^-J. 

Inserting this into (14.91) and using a > 1 and ( 14. 8ft 



4^5* + ^k 2 ) s 2 ds. 



o 



d 
d7 





2a 


"2w 2 /" 










. r Jo 
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47rjs a ds - Ew I I 47r5 + -k 2 s 2 ds 



( 



3Mw* 



M 1 



47rajs 3 ds 



(4.10) 



notice that the first term coming from the above estimate for a r is exactly 
cancelled by the last term in (14. 9p . which is crucial. Now the estimates 
w 2 /E < w and E > 1 imply that 



7T 



S(t,r)<j(t,r) + 



OO POO 



oo JO 



<j(t,r) + 



7T Li 



r 4 W. 
U 1 



E 

OO POO 



f(t, w, r, L) dL dw 



w f(t, w, r, L) dL dw 



W- r 2 



-oo JO 

j(t,r). 



(4.11) 
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By (12351) . re(t,r) = for r < i2_(t), and 



a+M 



47T /" r 4-7T f" 

\K(t,r)\<a + — j(t,s)s 2 ds < a + — p(t,s)s 2 ds < (4.12) 
r Jo r Jo 



Let 



J(t,r) 



/ 4vrj(t,s)s 2 rfs > 0. 
Jo 



Then using the fact that J is increasing in r and R-(t) > r_ for i e [0, io[> 



/cVds < a+M 



R-(t) 



4:WjT 2 dr I ds < 



a\M 
Rjt) 



J < 



a 2 M 



By the definition of 7, see (14. 5p . 

1 

r 



47rj(t, s) s 3 ds > jJ(t,r) 



J. 

(4.13) 
(4.14) 



We insert the estimates (ETTTj) . (fl~T3"l) . and (T4TT4|) into (flTTUD and use the fact 
that r > R-{t) > r_ along the characteristic under consideration to obtain 
the estimate 



d 
dt 



1 - 



2M 



> 2a j 



2jw 2 



1 + 



1 



W-ri 



Ew-- 



9a 2 ,M 



2 r_ 



3Mw z + — — 
r_ / r_ 



Due to the estimate Ew < E 2 = 1 + w 2 + L/r 2 this implies that 



dt 



'( 2M\ ~ 






K 1 - r) E \ 




[( 









9a 2 ,M 
27-I ± — 

' 2 r_ 



3a + M 



1 + 



+ 7C 



9 a 2 M 



2 r_ 



1 + 



If we recall the definition of 7 in ( 14. 5 j) we see that the condition ( I3.2p im- 
plies that the factor multiplying w 2 in the right hand side of this estimate is 
positive. Hence the estimate is still true if we replace w by W-, since by as- 
sumption w 2 > W 2 on the time interval [0, to[. Again recalling the condition 
( 13.21) we conclude that 



d_ 

dt 



1 - 



2M 



E 



> 
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as long as t G [0,to[- Hence for any characteristic in supp/ on this time 
interval, 

2M\ / „ s L \ ( 2M\ ( o . . L 
1 " "TIT 1 + w ( f ) + ~T77\ > 1 " -TTvT 1 + w (°) + 



r(t) y V ^ 2 0) / V r (°) / V r2 (°) 

Thus 

(*) + I?7TV > 1 - 177^ V) 2 (0) + 



r 2 (t) ~ V r (°) /V r2 (°) / r (°) ' 

Since f > 0, r(t) > r(0) > r_. Together with w(0) > this implies that 

9/ . / 2M\ „ 2M / L + 

= (i-0- 2 -^(i + ^)>^!, 

the latter because of ([S3]). Therefore w(t) > W_ for t G [0,t ]- From (Oil 
it follows that 

a+ f r . o , a+M 



Thus by ( 12321) . 



/0(*,r) > -— / A-Kjs 2 ds > 
r Jo 



a- W- a + M 

< r < 1 + — — . 

a+ ^/l + W 2 + L + /r 2 ~ ~ r_ 



Since r(0) e]r_,r + [, this implies that R-{t) < r(t) < R+(t) for t G [0,t ]; 
recall @3J and (|Q|) . Therefore gSJ holds on [0,t ], and t„ = T. 
Step 2: Global existence in maximal- areal coordinates. 

The next step in the proof of the theorem is to show that the estimates 
obtained in Step 1 on the maximal existence interval [0, T[ imply that T = oo. 

In [TUl Thm. 2.1] sufficient conditions for global existence in maximal- 
isotropic coordinates are obtained. It is straightforward to see that these 
conditions are also sufficient in our case, i.e., if we can show that the support 
of the momenta and the metric function a do not blow up, then T = oo fol- 
lows; the local well-posedness in maximal-areal coordinates and continuation 
of solutions will be studied more systematically in [3]. At points (t,r,9,<f>) 
with r < R_(t) there is no matter. Hence by fl2T5|) . ([2~20D . Q , and fl2TTTD 
this implies that 

a r (t,r) = a r (t,r) = n(t,r) = for r < R_(t). (4.15) 
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In particular the bounds in (14. Xp hold for all r > and t 6 [0, T[. In view of 
the latter it is sufficient to verify that 

sup{£ | (r,w,L) e supp/(t), t G [0,T[} < oo. (4.16) 

In what follows C > denotes a constant which only depends on the param- 
eters specified in ( 13. 2 j) and (13.31) and which may change from line to line. By 
( 12.71) and (I2.12D — (I2.14D it follows that along characteristics, 

• a r w an / 9 L \ . 

Since by ( 14. 6[) and (I4.3[) r(t) > R-(t) > C(l + t) for any characteristic in 
supp /, and g"2]) imply that 

C 

l«(*> 01 + <*r{t,r) < +t )2 - 

We point out that in this particular case when u> > the first term of the 
right hand side of ( 14.171) can be dropped since it is non-positive but this more 
general estimate will be needed later. Since a < 1 and a > 1 we thus obtain 
the estimate 

1*1 < C^, (4.18) 

and Gronwall's lemma implies the desired bound (14.161) . Thus T = oo. 

Step 3: Future geodesic completeness. 

In order to analyze this question polar coordinates such as we use on the spa- 
tial t = const surfaces are not convenient, and we introduce new coordinates 

x = t, x 1 = r sin 6 cos 0, x 2 = r sin 9 sin 0, x 3 = r cos 6, 

which can be thought of as the corresponding Cartesian coordinates. In these 
coordinates the metric becomes 

™ . nr* . nr* , 

S'oo = -« +a fj , g 0i = a p—, g {j = + {a - l)—^-, 

where Latin indices i,j run from 1 to 3 and X\ = bijxK Let us now con- 
sider an arbitrary future directed, time- like or null geodesic, i.e., a solution 
(x 7 (s),p 7 (s)) of the geodesic equations 

dx~< 

-P 1 , 



ds 
dp 1 

ds 
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where Greek indices 7, 5, e run from to 3, Tj e are the Christoffel symbols, 
and 



m 2 is conserved along any geodesic, and p° cannot change sign. Such a 
geodesic exists on a maximally extended interval [0, s + [, and future geodesic 
completeness means that s + = 00 for all such geodesies. 

The following relations between the variables r, w, L, and p 1 hold: 



That we had m — 1 before means that our system describes an ensemble of 
particles of rest mass 1, but for geodesic completeness we have to consider 
any time-like or null geodesic, i.e., we have to allow any m > 0. Since 
dt/ds = p° > we can re-parameterize the geodesic by coordinate time 
t E [0,t+[. Along the geodesic again fTCT7|) holds. By (ETT5|) . E = as 
long as r < R-(t) so that also (I4.18P holds along the geodesic, and E < C 
for some constant C. Since by (14. ip a is bounded from below, p° remains 
bounded on [0, t + [. Since both a and (3 are bounded as well also the p l remain 
bounded, and hence t + = 00. But since dt/ds = p° is bounded, this implies 
that s + = 00, and the proof of the theorem is complete. □ 
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